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We formulate a framework to determine the mass of glueball states of Landau gauge Yang-Mills 
theory in the continuum. To this end we derive a Bethe-Salpeter equation for two gluon bound states 
including the effects of Faddeev-Popov ghosts. We construct a suitable approximation scheme such 
that the interactions in the bound state equation match a corresponding successful approximation 
of the Dyson-Schwinger equations for the Landau gauge ghost and gluon propagators. Based upon 
a recently obtained solution for the propagators in the complex momentum plane we obtain results 
for the mass of the 0 ++ and 0 1 glueballs. In the scalar channel we find a mass value in agreement 
with lattice gauge theory. 


I. INTRODUCTION 

The ’physical’ spectrum of pure Yang-Mills theory 
is made out of glueballs [![. There is substantial ev¬ 
idence, both from lattice calculations (see S i and 
Refs, therein) and results from Dyson-Schwinger equa¬ 
tions (DSEs) for the Landau gauge gluon propagator [J- 
@ that transverse gluons violate positivity and therefore 
cannot be part of the asymptotic state space of the the¬ 
ory. Consequently, the first physical excitation of the 
Yang-Mills vacuum is the lowest lying glueball state. 

It is an important task to determine the mass of this 
state. Indeed, the properties of glueballs have been inves¬ 
tigated since their prediction in the middle of the 1970s 
[lj. Today, the glueball masses in pure Yang-Mills theory 
are known rather accurately owing to high statistics lat¬ 
tice calculations EH3- Unquenched lattice calculations 
are also available, although there are considerable un¬ 
certainties in the determination of unquenched glueball 
masses EMI- This is mainly due to severe problems 
with the signal to noise ratio, thus requiring large statis¬ 
tics. In principle, it is also not easy to disentangle states 
with large glueball components from states dominated by 
other constituents such as quark-antiquark pairs. Natu¬ 
rally, this problem has a counterpart in the experiments: 
a glueball cannot be distinguished from a meson by quan¬ 
tum numbers and masses only. The determination of 
the decay channels of a given state is therefore vital for 
its identification. Ongoing and new experiments such as 
BES III [15] and PANDA (TgI | have dedicated parts of 
their programs to the identification of heavy glueballs in 
the charmonium region and beyond. 

Alternative theoretical frameworks such as Hamilto¬ 
nian many body (l7l - [20| and strong coupling methods 
[2l| , potential approaches [2)1], Wilson loop based cal¬ 
culations [23|, flux tube models [23, chiral Lagrangians 
[26l-[28j], light front quantization [29| and the AdS/QCD 
approach f30j have shed some light on potential mass 
patterns and identifications of experimental states dom¬ 
inated by their glueball content. However, it seems fair 
to state that a detailed understanding of glueball forma¬ 
tion from the underlying dynamics of Yang-Mills theory 


is still missing. In this paper we report on further steps 
towards such an understanding. 

Working in Landau gauge, we construct homogeneous 
Bethe-Salpeter equations (BSEs) for glueballs which take 
into account the dynamics of gluon and ghost propaga¬ 
tion as well as their interactions with one another. This 
is detailed in Sections El and IIII1 where we also discuss 
the general form of Bethe-Salpeter vertices for any quan¬ 
tum number. The calculation is performed in Euclidean 
space, which implies that the bound-state consitutents 
are probed for complex momenta. In the literature, ex¬ 
ploratory BSE calculations using instantaneous approxi¬ 
mations mi or extrapolations of the propagators into the 
complex momentum plane [32j can be found. In this work 
we present first results for self-consistent and covariant 
BSE calculation of glueballs in Landau gauge using ex¬ 
plicit solutions for the ghost and gluon propagator DSEs 
in the complex momentum plane [6]. These are summa¬ 
rized in Section HVl 

II. BOUND STATE EQUATIONS FOR 
GLUEBALLS 

Our goal in this section is to provide a Bethe-Salpeter 
equation describing a glueball made from two gluons that 
are solutions of the DSE for the gluon propagator. A sim¬ 
ilar concept has proven very successful in the context of 
mesons, where the BSE of a quark-antiquark pair is used 
in connection with the corresponding DSE for the quark 
propagator, see e.g. [33l l35i| for reviews. A key property 
of this framework is consistency of the approximations 
made in the DSE and BSE. For mesons this implies to 
satisfy an axial Ward-Takahashi identity thus implement¬ 
ing constraints due to chiral symmetry and its breaking. 
One way of devising such a truncation is to derive both 
the truncation of the DSE and the truncation of the BSE 
on common grounds using a two particle irreducible ef¬ 
fective action (2PIEA). In the following we work along 
this strategy. Since we are working in Landau gauge 
we need to take into account the Faddeev-Popov ghosts. 
Thus we need a generalization of the usual BSE scheme 
that allows for mixing of bound states of different fields. 
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In the following we will give a derivation of a suitable 
set of bound state equations that provide the necessary 


couplings of bound state amplitudes with different field 
content (3(il|. 

We consider the following 2PIEA, 


P[D,G] = ^Tr In D 0 D~ l + ^TiD^D - Trln G 0 G _1 — TtGq 1 G + r 2 [Z), G], 


(1) 


where D and G are the gluon and ghost propagators. 
The interaction term is given diagrammatically by 


r 2 [AG] 



Each term contains one bare and one dressed vertex, 
the latter being represented by the shaded circles. The 
2PIEA is already truncated, i.e. we have left out all dia¬ 
grams including the four-gluon interaction. Furthermore, 
the dressed ghost-gluon and three-gluon vertices are as¬ 
sumed to be represented by suitable explicit expressions 


that capture the essence of the non-perturbative interac¬ 
tions. Such Ansaetze have been employed successfully in 
the past 0|; we come back to this point in Section IIVI 
The corresponding Dyson-Schwinger equations for the 
ghost- and gluon propagators can be found by variation 
of the effective action with respect to a propagator, i.e. 

5T[ ^ G] = ~ D~ l + D- 1 + S d [D, G } = 0 (3) 
6T f c G] = 2G" 1 - 2Gq 1 + E G [D, G] = 0, (4) 

where we have T,a = yf with A £ {D,G}. Diagram¬ 
matically, the resulting DSEs read 



We now proceed along the lines of Ref. fH]. In the fol¬ 
lowing we will use a shorthand notation omitting the 
space-time arguments and indicating primed arguments 
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by primed functions. We denote the solutions of DSEs 
© and fl]) by D and G and perform a variation in two 
variables. Keeping only the linear terms we arrive at 


5T[D, G] 

6D 

5T[D,G\ 

SG 


D-\-6d ,G-\-6g 


D-\-8 d ,G-\-6g 


ST[D,G\ 
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D,G 


<sr[£>, G] 
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d^x'cfiy' 


S 2 Y[D , G] 
5D5D' 



S 2 T[D , G] 
5GSD' 


5' d + 
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S' D + 
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S 2 T[D,G] 

5D8G' 


/ jA 
cat 


S 2 T[D,G] 

6GSG' 


S'n 


D,G 


S'n- 


D,G 


( 7 ) 

( 8 ) 


Using again the equations of motion we require for the 


solutions D and G to be stable that 

S 2 T[D,G] 


A ,a,S 2 T[D,G\ 


drx dry 


S', 
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D,G 


5DSD' 


SD5G' 


S' G = 0 

D,G 

( 9 ) 
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and 


/ 


d i x'd A y 


SGSD' 



S 2 T[D , G] 
5G5G' 


S' G = 0 . 

D,G 

( 10 ) 


Diagrammatically, the scattering kernels of the BSEs can 
be obtained by cutting a further line in the self energy 
diagrams with respect to the desired second constituent 
in the bound state. The variations 5' D ,5' G are identified 
with the Bethe-Salpeter vertices \d and xg- 

We then find the following coupled system of BSEs for 
ghost and gluon bound states 




-2 




where the arrow indicates symmetrization of the kernels 
with respect to the dressed vertices. The resulting cou¬ 
pled system of two-body equations serves to describe 
glueballs as bound states of either a gluon or a ghost- 
antighost pair. The latter is necessary in Landau gauge 
and represents contributions from the Faddeev-Popov de¬ 
terminant to the glueball masses. We will later discuss 
the relative importance of both contributions in differ¬ 
ent channels. For now we just emphasize that neither 
the ghosts nor the gluons are physical constituents in the 
sense that they do not appear as propagating particles in 
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the positive definite part of the asymptotic state space 
of QCD M. Note, that there is no mixed gluon-ghost 
(or gluon-antighost) contribution to the glueball vertex. 
Such gluon-(anti)ghost bound states, if existent, may be 
members of a BRST quartet together with transverse glu- 
ons and, thus, part of the unphysical Hilbert space [37j. 

The above system of BSEs within Yang-Mills theory 
can be further generalized to full QCD by including quark 
contributions. Considering the corresponding effective 
action 


r[D,G,S] =^Trln D 0 D _1 + ^Tr D^D 

— Tr In GoG" 1 -TrG^G 

— Tr In SqS -1 — Tr Sq 1 S + T 2 [D, G, 5], 

(13) 


with diagrammatic representation 


r 2 [AG,s] 



we can apply the same derivation as before considering 
variations with respect to all types of propagators. We 
then find the full system of coupled bound state two-body 
equations 



first of (flTl) provide for glueball/meson mixing. Although 
in this work we will restrict our explicit calculations to 
pure Yang-Mils theory we would like to add some com- 


This set of BSEs describes mesons and glueballs in an 
approximation that can be seen as a generalized ladder 
truncation. Note, that the last diagram in (1151) and the 
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ments on the influence of these terms onto the flavor sin¬ 
glet meson spectra. 

In the pseudoscalar channel these terms generate a con¬ 
tribution to 77 - 77 ' splitting. In the framework of BSEs 
the conventional approach to this problem is to include 
beyond rainbow-ladder terms connected with the axial 
anomaly in the quark-gluon interaction (38l - f40l| . Our 
framework provides the additional effect of a direct mix¬ 
ing of the flavor-singlet r^-meson with the J PC = 0 *" 
pseudoscalar glueball. Both effects together affect the 
? 7 o-meson mass. Our set of BSEs provides the simplest 
means of consistently including such glueball/meson mix¬ 
ing into BSE calculations. However, since in lattice cal¬ 
culations the mass of the pure glue pseudoscalar glueball 
is found to be around 2.5 GeV [9], it is not clear how large 
these mixing effects might be. 

Considerable mixing effects may be expected, however, 
in the scalar meson sector, where there seem to be more 
states than one can accommodate in conventional quark 
anti-quark multiplets (JO]. While some of these states 
can well be accounted for by large four-quark components 
(42l - |44| , others may very well be characterized by a dom¬ 
inant glueball contribution. The set of Equations (fT51) to 
m may well provide a viable starting point for sophisti¬ 
cated investigations of a realistic scalar meson spectrum. 

In the following, we will focus on the coupled system of 
bound states for a pure gauge theory neglecting quarks 
(Eqs. (1111) and (fT!?l) l. To solve this system numerically, we 
need reliable information on the nonperturbative propa¬ 
gators of ghosts and gluons as well as a general expres¬ 
sion for the bound state vertices \d and \G- In the next 
section we will discuss the latter, providing suitable ex¬ 
pressions for arbitrary quantum numbers of glueballs. 


III. BOUND STATE VERTICES FOR 
GLUEBALLS 

We will now show how suitable bound state vertices 
Xd and xg to be used in Eqs. (ED and m can be con¬ 
structed. We start our present discussion from some gen¬ 
eral observations. 

A bound state of two relativistic particles can be de¬ 
scribed by three quantum numbers: total spin J, parity 
P and charge parity C. Furthermore there are only two 
characteristic momenta involved, conveniently chosen to 
be the total momentum t^ and the relative momentum 
r^. These two vectors can be used to construct suitable 
vertices xd and \G for our bound state problem. The 
idea is to construct basic invariant vertices with correct 
parity and charge conjugation properties and supplement 
these with an appropriate tensor representing a given to¬ 
tal spin. 

Let us first consider the vertex xd for the two-gluon 
bound state. A general scalar bound state vertex has to 
transform like a rank two Lorentz tensor for gluonic con¬ 
stituents, yet it has to be invariant under Lorentz trans¬ 
formations for it is meant to represent a scalar. Thus we 


are looking for a structure that transforms like 

A%A \T kX = T^. (18) 


There exist only two tensor structures with J = 0 that 
satisfy (1181) . namely the metric tensor and a combina¬ 
tion of the totally antisymmetric tensor and characteris¬ 
tic momenta: 

= r ^tx (19) 


Whereas the first represents a parity-even state, the sec¬ 
ond choice is odd under parity transform. 

For the vertex xg composed of ghost fields the situa¬ 
tion is trivial, since we are looking for a term that cou¬ 
ples to scalars and transforms like a scalar itself. The 
appropriate vertex is the identity in Lorentz space and 
has positive parity. This restricted choice has interesting 
implications as discussed below. 

In addition to the basic tensors in Eqs. m, represent¬ 
ing the Lorentz structure of the constituents, we need 
suitable tensors representing a given total momentum J 
of the bound state. A Lorentz tensor representing a mas¬ 
sive field with total spin J is required to have precisely 
2 J 4- 1 independent components to represent the possi¬ 
ble spin polarizations. The construction of such tensors is 
known and a detailed treatment can be found e.g. in (46j . 
We will repeat parts of the construction here in a slightly 
more explicit form focused directly on the construction of 
Bethe-Salpeter vertices. Consider first tensors T ai ,...,aj in 
three-space of rank J. To represent angular momentum 
J we require the tensor to be symmetric in all indices 
T ai ...aj = T'p[ ai ...aj] and traceless with respect to any 
pair of indices ]T) m = 0. The first constraint 

leaves the tensor with i (J 1 2 + 3 J + 2) independent com¬ 
ponents, while the second one imposes ^ ( J 2 — J ) further 
restrictions, thus leading to a tensor with 2J + 1 inde¬ 
pendent components, as required^ The construction of 
tensors in three-space is now easily transferred to four- 
tensors. If we require the tensor to be transverse 

to the total momentum of the particle in every index 
t v T... v ... = 0 and adopt the particles rest-frame, we see 
that all components with time-like indices vanish, leaving 
only components with space-like indices. So we are left 
with nothing else but the three-tensor considered before, 
which has 2 J + 1 independent components. Thus we find 
the constraints for a Lorentz-tensor T,/ ,,. of rank J to 

represent angular momentum J: 


1. T is symmetric in all indices, 


rjij 

1 


= T. 




( 20 ) 


1 It is also possible to construct tensors representing half-odd in¬ 

teger spin. Since we are dealing with Bethe-Salpeter equations 
of two particles in the same representation, so that the total an¬ 
gular momentum is integer, we will not consider this possibility 
here but instead refer the interested reader again to 14611 . 
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2. T is transverse to the total momentum of the par¬ 
ticle in every index, 

t v T J V = 0. (21) 


3. T is traceless in every pair of indices in the rest- 
frame, 


T" / ; a ; a - = 0. (22) 

For the glueball masses we only need one such tensor from 
each multiplet. To construct such a tensor for angular 
momentum J one can build the J-fold tensor product 
of a transverse projector that transforms like a vector 
and then subtract the traces with respect to every pair 
of indices. Starting with J = 1, a suitable transverse 
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four-vector can be obtained by contracting the transverse 
projector t™ (with respect to the total momentum t ) and 
the relative momentum r, 


Qu 



VM r v 
t 2 3 J 


_ (r-t) tA 

t 2 )' 

(23) 


With only one Lorentz index this transverse vector al¬ 
ready gives a possible angular momentum tensor for 
J = 1. For higher J one builds symmetric J-fold ten¬ 
sor products of (E51) by 

= Q mi * Qfij- ( 24 ) 

The next step is to remove the traces of these tensors with 
respect to every pair of indices. This can be achieved with 
the general formula 


T u 


“ (2J- 1) _1 Y, 


Q t 




where denotes the sum over all essentially different 
permutations of the indices @ We furthermore define 


h=T 2 - 


( r ' tf 

t 2 


and the tensort@ 

f>J,j _ p r r /) (J — 2j) 


{M1M2 /xj} ) 


(26) 


{/2l/22"-°A‘2 3 -l ^ 2j < J, 

(27) 

R J J = f J / 2 X, 


2 j = J. 
(28) 


Using E5l) and EH) together with the above definition, 
we finally obtain the desired total spin tensors in closed 
form as 


...pj — 

2j< J / j \ ^ 

e (n 2^ - /.-)+i) 

j =i \fe=i / 

(29) 

With these we have access to higher orbital angular mo¬ 
mentum states built for glueballs with two gluon con¬ 
stituents (with C = +1) in the Lorentz singlet channels. 


2 This means that the sum has to be divided by appropriate com¬ 
binatorial factors. 

3 For the convenience of the reader we have denoted the rank of 
the raw tensors as a superscript. 


(2 J 1) (2 J 3) r MlM2 T M3M4 Q hzXfls [Aj 

p 
1 P- 


(25) 


For arbitrary even J we can use, 
r £^ 1 ...»At*y,0)= T ^ J A(t 2 ,r 2 ,0)g ltv , (30) 

r )Ljn... M .,(*V 2 ,0) = T m-HJ A (t 2 ,r 2 ,e) (r-t) r K t x e KX ^, 

(31) 

and if J is odd, 

r £^i = T lll ... llJ A(t 2 ,r 2 1 8) (r-t) r K t x e KXliU , 

(32) 

=T fll ...» J A(t 2 ,r 2 ,0)g tlv . (33) 

Here we have introduced scalar functions A(t 2 ,r 2 ,8), 
which are even under inversion of the angle 8 between 
r and t. The additional factors (r • t) ensure the correct 
behavior of the bound state vertices under charge parity 
transformations, which result in a simple flip of the sign 
of the relative momentum r^ in our framework. 

The corresponding vertices for glueballs with a ghost 
anti-ghost pair as constituents are constructed along the 
same lines. For arbitrary positive parity even J we can 
use 

CL,(*V 2 -0) = T ^j A (t 2 ,r 2 ,e), (34) 

and if J is odd the negative parity states are obtained 
from 

r iZ^y^)= T ^-^A(t 2 y,8). (35) 

Note that the parity P of the angular momentum tensors 
Tm-.-ijj is given by P = (—\) J . Thus the parity of the 









6 


even (odd) J tensors is positive (negative). Hence there 
are neither contributions from ghost anti-ghost pairs to 
glueballs for even J and with quantum numbers J nor 
for odd J with J ++ . For the gluonic vertices the natural 
parity of the total spin tensors can be supplemented by 
the odd parity basis element given in (fl9ll , thus changing 
the overall parity of the tensor representation. This is 
not possible for the ghost vertices. Consequently, we only 
find the restricted set ( 1551 ) and ( 1551 ) of possible quantum 
numbers for glueball states with ghost contributions. 

Another potential restriction for the contribution of 
gluonic vertices to glueballs has been frequently discussed 
in the literature in the context of model building [47| : 
if the gluonic constituents were massless and on-shell, 
Yang’s theorem (48) would restrict the number of al¬ 
lowed quantum numbers drastically. In our framework 
this constraint appears to be almost irrelevant. The non- 
perturbative gluonic constituents that appear in the BSE 
are neither on-shell nor massless. Instead they acquire a 


dynamically generated mass, as discussed in more de¬ 
tail in the next section^ Thus, gluonic contributions 
in Yang-forbidden channels may be suppressed but cer¬ 
tainly not forbidden. These channels are the 1 ++ and all 
odd J ^-channels. Since in the l ++ -channel ghost con¬ 
tributions are absent as well, the potentially suppressed 
gluonic contributions may lead to an ‘unnaturally’ large 
glueball mass in this channel. This is indeed observed 
in lattice calculations [g, 9j|. In the odd J ^-channels, 
however, ghost contributions are allowed. If the suppres¬ 
sion of the gluonic contributions were strong, these states 
could be termed ‘ghostballs’. We will study such states 
in future work. 

In addition, there may be another basic vertices for 
the gluons in the spin 2 channel, traceless symmetric ten¬ 
sors. For example, Landau constructed one for QED with 
J PC = 2 ++ [45]. The explicit form of this basic J = 2 
tensor with our notations would read 


F 2 = t 4 l- 


/I i 1 \ 

1 /J 2 + ^9ihh9h^2 + 9 ii ii2 9 v hi J 

2 (l 1 1 1 1 1 \ 

t y—ginitfii^n 2 — 7^9nni^vtH2 — — — 9 9v A J 


1 M2 > 


(36) 


where the pair (p, v) denotes the Lorentz indices of the 
gluon constituents and (/xi,/^) the one of the bound 
state. With such a tensor vertices for a given set of quan¬ 
tum numbers J PC can be constructed in a similar way 
as from the two singlet tensors in Eqs. 

Having discussed the derivation of the bound state 
equations we will use to describe glueballs and the form 
of the necessary bound state vertices, we will now turn 
to the numerical part of our investigation. We will find 
that it is necessary to solve the system of coupled DSEs 
of ghost and gluon fields for complex momenta in or¬ 
der to use the resulting propagators in our calculation of 
glueballs. 


4 Note that in the context of this discussion of Yang’s theorem 
it is irrelevant whether the mass generation mechanism leads 
to an infrared vanishing gluon propagator (‘scaling’) or an in¬ 
frared finite propagator (‘decoupling’). Furthermore, ’mass gen¬ 
eration’ in this context does not mean that the gluon propaga¬ 
tor acquires a pole at time-like momenta, but merely that the 
zero-momentum pole of the free propagator disappears due to 
interactions. 


IV. THE YANG-MILLS SYSTEM IN THE 
COMPLEX PLANE 


We solve the coupled system of bound state equations, 
m and (fT51) in pure Yang-Mills theory in Euclidean 
momentum space. In the rest frame of the glueball, 
its total momentum is then given by (0, 0,0, ims) with 
ms is the bound-state mass. Without loss of general¬ 
ity, the total momenta can be shared equally between 
the two constituents. Their momenta are then given by 
r± = (r ± f) /2, with relative momentum r between the 
constituents. It is then clear that the internal propaga¬ 
tor lines in BSEs are given by solutions of the DSEs for 
complex momenta. 

These can be obtained from the corresponding cou¬ 
pled set of Dyson-Schwinger equations. In Landau gauge 
the ghost propagator Dc{p 2 ) and the gluon propagator 
D^ip 2 ) are given by 


D^{P 2 ) 

d g (p 2 ) 




G( P 2 


PhPu\ Z(p 2 ) 


(37) 

(38) 


where the diagonal color structure has been omitted for 
brevity. Note that the gluon is transverse also non- 
perturbatively, thus spurious glueball states due to longi- 










Figure 1. Left: Results for the gluon dressing function Z(p 2 ) from the DSEs [gl] for real and space-like momenta, compared 
with lattice calculations ffiTI | . Right: Results for the gluon spectral function from DSEs for time-like momenta. Shown is the 
result from Ref. Q together with result obtained in the same truncation scheme but with slightly different numerics, see main 
text for further explanations. 



Figure 2. Real part of the gluon propagator function D(p 2 ) 


Z(p 2 )/p 2 in the complex z = p 2 plane jtf]. 


tudinal modes, as present in some potential models [22]. 
are naturally avoided. The coupled system of DSE (omit¬ 
ting two-loop diagrams) has been displayed diagrammat- 
ically in Eqs.®)-®. This system of equations has been 
considered frequently in the past years. It has been 
solved analytically in the deep infrared, where exact solu¬ 
tions without any truncations are possible [49|, [50]. Two 
qualitatively different solutions have been found named 
‘scaling’ and ‘decoupling’ (sEIHI]- Whereas the scaling 
solution consists of infrared power laws for all Green’s 
functions with an infrared vanishing gluon propagator 
and an infrared divergent ghost, the decoupling solution 
is characterized by an infrared finite gluon propagator 
and a finite ghost dressing function. Current lattice cal¬ 
culations on large volumes clearly favor the decoupling 
type of solutions (53]; there is, however, an ongoing dis¬ 
cussion on potentially significant effects from different 
gauge fixing strategies in the deep infrared (UE3- In 
this work we concentrate on decoupling type of solutions, 
which have been associated with a dynamically generated 
‘gluon mass’ll As already mentioned above, these still 


5 Of course, this ‘mass’ is not to be identified with the mass of a 


maintain transversality. At finite momenta, the equa¬ 
tions have to be solved numerically and approximations 
for the dressed vertices need to be introduced. A suit¬ 
able truncation scheme has been introduced in Ref. |60] 
and improved in Ref. [51] . It involves educated ansaetze 
for the ghost-gluon and three-gluon vertex and neglects 
the effects of the four-gluon interaction completely. The 
resulting modified system of equations has exactly the 
structure of the system ® and a obtained from the 
variation of the 2PIEA m in Section El For real mo¬ 
menta, we show the numerical solution for the gluon 
dressing function in Fig. [TJ As can be seen from the 
comparison with the lattice results [61], there is very 
good agreement in the infrared and ultraviolet momen¬ 
tum region, whereas in the mid-momentum region one 
observes quantitative deviations. These deviations are 
certainly in parts due to the neglected four-gluon inter¬ 
actions in the DSEs ®-®. They can be compensated, 
however, by suitably optimizing the input used for the 
dressed three-gluon vertex by simultaneously solving its 
own DSE together with those for the propagators [62j ■ 


physical particle. The analytic structure of the gluon propagator 
is clearly different from a simple mass pole [(J. 


























The system of DSEs (0, © for the ghost and gluon 
propagator has been solved in the complex p 2 plane re¬ 
cently, see Ref. jfii] for details. As explained above, this 
complex solution constitutes a vital input into the corre¬ 
sponding Bethe-Salpeter equation for the glueballs and 
is used in the following. In this respect it is important to 
note that the analytic structure of the gluon and ghost 
propagators as obtained in Ref. [6] shows branch cuts 
along the time-like momentum axis, i.e. for negative in¬ 
variant momentum squared, but no singularities away 
from the real axis in the complex p 2 plane. For the real 
part of the gluon propagator this can be seen in Fig. [5] 
corresponding plots for the imaginary part and the ghost 
dressing function can be found in Ref. Q. This behav¬ 
ior greatly helps in the numerical treatment of the BSE. 
Along the time-like axis of negative p 2 one can extract 
the gluon spectral function, which is shown in the right 
plot of Fig. [T] We have plotted the result from Ref. j(| 
together with a corresponding result obtained from an 
improved numerical method. The corresponding results 
are very similar except on a narrow region around the cut 
on the negative p 2 -axis. As a result, one obtains a con¬ 
siderable smoother spectral function, as can be seen in 
Fig. m We use this improved result in the present work. 

Finally, a comment on the scale is in order. This is 
fixed by comparison with the lattice results for the gluon 
propagator and remains fixed, i.e. it is not adjusted again 
in the glueball calculations. Thus, in principle, we obtain 
absolute values for the glueball masses. 


V. LOWEST LYING GLUEBALL MASSES 

We have solved the BSEs for a glueball in the scalar 
and pseudoscalar channel using the bound state equa¬ 
tions ED and ED together with the vertices (1501) to 
El- For the propagators of ghosts and gluons we use 
the numerical results discussed in the last section. 

Bound state vertices are not primitively divergent ver¬ 
tices and therefore they generically go to zero like power 
laws for large momenta [6.‘il |. In contrast to the behavior 
of meson BSEs, however, ED and ED also admit so¬ 
lutions with a logarithmic behavior in the UV. (Details 
on the asymptotic behavior of glueball BSEs together 
with an explicit analytic analysis will be given elsewhere.) 
These solutions do not correspond to bound states. In or¬ 
der to guide the iterative numerical procedure to the cor¬ 
rect bound state solution it turns out to be sufficient to 
introduce an additional Pauli-Villars term into the purely 
gluonic diagram of ED that depends on the momentum 
of the exchanged gluon propagator. Namely, we replace 
Z(k 2 ) —> Z(k 2 ) (l + k 2 / Ap V .) 1 , with k 2 the momentum 
of the exchanged gluon and A py a cutoff scale. By in¬ 
spection of the Bethe-Salpeter vertex functions we have 
verified, that such a term does not simply modify the log¬ 
arithmic solutions of the BSE above the scale Apy but 
indeed drives the equation to a different and well-behaved 
solution. We find that the resulting glueball masses are 



masses (GeV) 

J PC 

lattice 

Hamiltonian/ 
Regge theory 

constituent 

models 

this work 

0++ 

1.71 (5)(8) [9] 
1.73 (5)(8) [8] 

1.98 [18] 

1.58 [23] 

1.71 [25] 
1.86 [22] 

1.64 

(T+ 

2.56 (4)(1) [9] 
2.59 (4) (13) [8] 

2.22 |T8] 

2.56 ]23] 

2.61 [25] 
2.49 [22] 

4.53 


Table I. Scalar and pseudoscalar glueball masses (in GeV) 
from various studies. We quote the Model B data from 
Ref. gH. 

insensitive to all values A 2 PV > 100 GeV 2 of the scale 
that we have probed. 

We present our results in Table Q] together with corre¬ 
sponding ones from lattice gauge theory, the Hamiltonian 
approach and Regge theory. Additionally, we compare 
to a rather recent calculation in a non-relativistic con¬ 
stituent model. 

Comparing with the lattice results, we find that the 
state with quantum numbers 0 ++ is well reproduced on 
the five percent level. Compared with the lattice, the 
good agreement of our result for the lowest lying scalar 
glueball is remarkable, though probably not surpris¬ 
ing. As explained above, our truncation scheme for the 
ghost/gluon DSEs produces solutions which are point- 
wise similar to the lattice results in the low and high mo¬ 
mentum region and display a twenty percent difference 
for momenta around 1 GeV. Thus the overall quality of 
the truncation scheme is well below the twenty percent 
range and thus in agreement with our findings for the 
scalar glueball mass. The remaining deficiencies in our 
truncation scheme are in the details of the three-gluon 
and the missing four-gluon interactions. 

In contrast, the mass of the pseudoscalar glueball is 
much higher than that predicted by lattice calculations 
as well as by other approaches. As discussed in Section 
|m| there are no ghost contributions in these channels, 
leaving a greatly reduced BSE with only one gluonic di¬ 
agram to be solved. This diagram is, in turn, largely 
dominated by the three-gluon vertex, both directly and 
via the solution of the DSE for the gluon propagator. 
Since the ansatz used here was devised in the context 
of the study of gluon and ghost DSEs for real momenta, 
it is conceivable that the behavior of this ansatz in the 
complex plane affects significantly the glueball spectrum 
in particular for states with no ghost-antighost content. 
The study of the connection between the details of the 
non-perturbative gluon self-interactions and their impact 
on glueball masses will be the subject of future work. 

VI. SUMMARY 

In this paper we have presented a framework that al¬ 
lows to calculate glueball properties from the dynamics 
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of Landau gauge Yang-Mills theory. We have constructed 
a set of bound state equations that includes both ghosts 
and gluons degrees of freedom thus taking into account 
also the effects of the Faddeev-Popov determinant. This 
set of equations allows for mixing of bound state contri¬ 
butions from different species of particles and is readily 
generalized to full QCD, including quarks. It thus nat¬ 
urally incorporates meson/glueball mixing. Furthermore 
we have presented suitable representations for the bound 
state vertices for arbitrary quantum numbers J PC . 

As an illustration of the framework, we have calcu¬ 
lated the scalar and pseudoscalar glueball mass. Our 
result for the scalar glueball state is certainly encour¬ 
aging, although in the pseudoscalar channel the mass is 
exceedingly high. Compared to the recent exploratory 
approach of Ref. j32 | we have made a number of techni¬ 
cal improvements. Most important are the use of explicit 
solutions of the ghost and gluon propagators in the com¬ 
plex momentum plane. Furthermore, our approach fully 
maintains multiplicative renormalizability. 

Our framework for calculations of glueball properties 
offers various prospects of improvements and applications 


in the near future. First, different ansaetze for the three- 
gluon vertex should be used and its impact on the spec¬ 
trum analyzed. It would be desirable, although techni¬ 
cally very demanding, to use dynamical three-point ver¬ 
tices as in [62| and to include the four-gluon interaction 
contributions fodl . |65| into our framework. On the other 
hand, a very important extension is the inclusion of me¬ 
son/glueball mixing along the lines of Eqs. (fl5l) to (TT71) . 
This will allow to leave the sector of pure gauge field cal¬ 
culations of glueballs and thus provide access to realistic 
glueball properties in the future. 
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